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ABSTRACT 


A  regional  event  discriminant  is  being  developed  and  tested  based  on  the 
waveform  of  high  frequency  The  data  base  being  used  in  Che  development 
consists  of  signals  from  explosions  and  earthquakes  recorded  on  the 
western  U.S.  digital  network.  This  net  consists  of  the  four  LLNL  stations 
MN\',  KNB,  ELK  and  LAC;  the  two  DWWSSN  stations  ALQ  and  JAS  and  the  two 
university  run  stations  PFO  and  PAS.  'It  has  been  discovered  that  at 
most  stations  the  waveform  of  P^  onset  is  remarkably  stable  and  different 
from  the  corresponding  waveforms  from  earthquakes.  A  forward  modeling 
study  of  broad  band  explosion  P^'s  revealed  that  the  distinctive  character 
of  their  waveform  is  caused  by  a  strong  pP^  arrival.-  Depth  phases  from 
earthquakes  arrive  much  later  in  the  signal.  It  was  found  that  a  clear 
effective  pP^  arrival  was  present  in  all  cases.  However,  it  consistently 
arrives  later  than  the  predicted  elastic  time.  For  Pahute  events,  the 
amplitude  of  effective  pP  is  close  to  the  elastic  predictions.  For  Yucca 
Valley,  it  is  consistently  larger  Indicating  the  effect  of  a  site  dependent 
nonlinear  process  in  the  source  region.  An  appropriate  value  of  t*  for 
P,,  appears  to  be  in  the  range  of  0.1  to  0.2  seconds.  The  frequency 
content  of  the  explosion  P^^  energy  indicates  that  it  is  caused  by  turning 
ravs  in  the  lid  gradient  rather  than  true  head  waves  traveling  on  the 
crust  mantle  interface.  The  discrimination  capacity  of  the  waveform 
was  measured  quantitatively  bv  corrfilating  the  average  explosion  Retrace 
with  a  data  base  of  explosion  and  earthquake  signals.  The  populations 
separated  to  a  significant  level  down  to  magnitudes  less  than  4.0.  It 
was  found  that  the  average  explosion  waveform  from  one  station  could  be 
used  to  discriminate  data  from  a  different  station  establishing  that  the 
Pp  waveform  discriminant  is  transportable. 
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INTRODUCTION 


High  frequency  regional  data  appears  to  be  among  the  most  unstable  and 
difficult  to  characterize  in  seismology.  For  this  reason,  most  of  the  attempts 
to  solve  the  regional  discrimination  problem  have  been  based  on  highly  empirical 
approaches.  Ratios  of  phase  amplitudes  in  either  the  time  or  frequency  domain 
have  been  tested  to  find  whether  they  are  systematically  different  for  earthquakes 
and  explosions.  Recently,  Taylor  et  al.  (1988)  have  demonstrated  an  impressive 
capability  for  event  discrimination  using  spectral  ratios  of  Lg.  In  essence, 

;  ley  demonstrated  that  regional  signals  from  explosions  are  systematically  lower 
in  frequency  content  than  earthquakes.  This  confirms  an  earlier  result  of  Murphy 
and  Bennett  (1982) .  Though  development  of  this  discriminant  represents  a 
significant  advance,  there  remain  important  questions  of  the  transportability 
of  the  method.  The  reason  is  that  a  deterministic  model  for  the  change  of 
spectral  behavior  between  earthquake  and  explosion  sources  is  not  available. 
There  is  no  evidence  in  teleseismic  data  that  explosion  sources  have  a  different 
spectral  decay  rate  than  earthquakes.  Most  investigators  assume  that  both  sources 
fall  off  as  f'2,  though  Burger  et  al .  (1987)  have  shown  that  this  is  not  necessarily 
true.  Studies  which  have  indicated  that  explosions  are  actually  enriched  in 
high  frequency  with  respect  to  earthquakes  include  those  of  Savino  et  al .  (1980), 
Everendon  et  al .  (1986)  and  von  Seggern  and  Rivers  (1979).  Thus,  if  the  spectral 
discriminant  works  consistently  on  Lg  in  the  western  United  States,  it  must  be 
related  to  phenomena  associated  with  Lg  propagation  there.  We  can  not  rely  on 
the  .same  processes  to  occur  in  a  consistent  fashion  in  other  regions. 

To  avoid  such  difficulties  in  transporting  discriminants  to  new  tectonic 
and  .seismic  r^c^imos,  it  will  be  necessary  to  develop  discriminants  based  on  a 
sound  understanding  of  the  physical  processes  which  underlie  them.  This  knowledge 
will  permit  an  assessment  of  whether  and  how  the  discriminant  should  be  modified 
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for  other  areas.  The  investigations  discussed  here  are  directed  at  establishing 
what  the  most  stable  portions  of  the  regional  P  wave  signal  are,  and  if  they 
can  be  modeled  in  a  deterministic  fashion.  Specifically,  we  will  discuss 
ohi;ervations  of  NTS  nuclear  events  at  stations  in  the  western  U.S.  digital 
network.  We  have  used  arrays  of  events  to  form  seismic  seccions  at  each  station 
and  examined  them  for  stable  arrivals.  Most  of  the  P„  and  Pg  wave  trains  exhibit 
a  low  level  of  coherence.  However,  we  have  found  that  the  very  onset  of  P,,  is 
remarkably  stable  at  high  frequency  and  that  it  has  a  character  strongly  influenced 
by  free  surface  phases.  Since  free  surface  phases  occur  later  for  earthquakes 
than  explosions,  it  is  possible  to  develop  a  discriminant  based  on  onset. 
Because  the  physics  behind  the  discriminant  is  straightforward,  transporting  the 
technique  to  other  regions  should  not  pose  a  significant  problem. 

For  the  purposes  of  this  discussion,  we  will  define  the  western  U.  S.  digital 
array  as  including  the  four  LLNL  stations  MNV,  KNB,  LAC  and  ELK;  the  two  DWWSSN 
stations  ALQ  and  JAS  and  the  Caltech  and  UCSD  stations,  PAS  and  PFO.  In  the 
following  report,  we  will  discuss  the  seismic  sections  of  the  regional  P  wave 
train  we  have  assembled  for  each  station  to  date.  We  will  compare  the  relative 
coherency  of  Pj.  and  Pg,  and  determine  the  average  P„  onset  waveform  by  stacking. 
We  will  identify  stations  with  complex  receiver  functions  and  those  being  affected 
by  complex  structure  in  the  lithosphere.  This  structure  explains  some  of  the 
in.stability  of  the  ratio  of  the  amplitude  of  the  P^  and  Pg  arrivals.  We  will 
then  present  modeling  studies  of  the  observed  P^  onset  waveforms  which  delineate 
the  effects  of  pP.  We  will  also  model  some  of  the  clear  effects  of  source 
scaling.  Finally,  we  will  demonstrate  how  P^  onset  can  be  used  to  discriminate 
explosions  from  earthquakes. 
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THE  UESTERN  U.S.  DIGITAL  ARRAY 

The  shore  period  instruments  at  the  two  DWWSSN  stations  are  relatively  band 
limited  in  that  they  are  designed  to  mimic  the  short  period  Benioffs  of  the  old 
UVSSN.  The  rest  of  the  stations  have  modern  broad  band  instruments.  To  make 
the  following  comparisons  uniform,  we  have  converted  all  of  the  data  to  the  short 
period  WVSSN  response.  In  later  sections,  we  will  consider  deconvolved  broad 
band  data.  Figure  1  is  a  map  displaying  the  stations  of  the  digital  array  as 
they  distribute  about  NTS.  Azimuthal  coverage  is  relatively  good,  and  stations 
Lac  and  PFO  are  almost  exactly  aligned  with  NTS.  Shown  as  circles  are  the 
locations  of  earthquakes  we  have  used  for  discrimination  studies.  The  closest 
station  is  MlFv'  at  a  range  of  195  km  to  the  nearest  NTS  events  and  the  farthest 
is  .AI15  at  890  km.  The  other  stations  are  at  about  325  to  375  km  in  range  on 
average . 

KN'V  The  first  station  in  the  array  we  wish  to  discuss  is  Mina,  Nevada.  As  shown 
in  Figure  1,  it  is  located  northwest  of  the  test  site.  Figure  2  shows  the 
vertical  component  observations  of  4  Pahute  events  at  that  station.  In  this  and 
all  of  the  following,  seismic  sections  are  shown  reduced  to  the  apparent  velocity 
of  P„;  that  is  the  arrival  time  is  selected  manually  and  a  ten  second  leader  is 
added  to  it.  The  P^,  barely  emerges  from  the  front  of  the  Pg ,  but  it  has  a  stable 
form  which  we  will  show  in  the  following  to  be  very  significant.  This  shape  can 
be  seen  most  clearly  in  the  bottom  trace.  The  arrival  has  a  small  upswing 
first  which  we  will  refer  to  as  the  "a"  swing.  This  is  followed  by  a  strong  "b" 
swing  downward  and  a  structured  "c"  swing  upward.  The  "c"  swing  exhibits  a 
double  shoulder  or  splitting.  In  the  following,  we  will  show  that  this  feature 
is  causf'd  by  the  arrival  of  pP,,. 
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Figure  1.  Tiie  western  United  States  digital  array  as  defined  for  the  purposes 
of  this  study.  MNV,  KNB,  LAC  and  ELK  are  LUIL  stations.  ALQ  and  JAS  are 
DWUSSN  stations  and  PAS  and  PFO  are  university  run  stations.  The  open  circles 
are  the  locations  of  earthquakes  In  the  discrimination  data  base. 


RANGE  (KM) 


Figure  2.  A  seismic  section  of  four  Pahute  explos 
is  the  small  signal  consisting  of  only  3  or  4 
larger  P  The  P^^  waveform  Is  stable  whereas  no  pori 


The  Pg  arrival  Is  as  unstable  ir.  character  as  most  high  frequency  regional 
signals  usually  are.  There  may  be  some  periodic  sequence  of  pulses  contributing 
to  the  upper  traces  which  could  be  caused  by  crustal  resonance  phases.  In  a 
previous  report,  we  associated  such  phases  with  successive  reverberations  of 
energy  in  the  cri:stal  wave  guide  (Burdick  at  al . ,  1988).  At  any  rate,  it  is 
clear  that  there  is  no  feature  in  the  regional  P  wave  signal  as  stable  as  the 
shape  of  the  P  wave  onset.  It  is  important  to  note  here  that  the  relative 
amplitude  of  to  P^,  is  about  5  to  1  for  the  purposes  of  the  following  discussion. 

Figure  3  shows  a  comparable  section  of  signals  from  Yucca  Valley  events. 
The  Pafiute  events  are  at  an  average  of  aoout  200  km  in  range,  and  the  Yucca 
events  are  at  about  231  km.  Thus  we  see  some  significant  P^  move  out  from  Pg. 
The  splitting  of  the  c  swing  is  apparent  in  all  traces,  particularly  the  bottom 
ones.  The  top  two  traces  show  a  very  deep  split,  but  given  that  this  is  high 
frequency  data  the  shape  of  onset  is  very  stable.  The  Pg  shows  no  stable 
features  in  this  format.  No  separate  phases  can  be  discerned  moving  through  the 
section  except  for  a  crustal  S  w^ve  at  the  very  back.  It  is  interesting  to  note 
that  much  more  S  is  generated  at  Yucca  Valley  than  at  Pahute  Mesa.  Presumably, 
this  is  ciue  to  the  stronger  scattering  at  Yucca  Valley. 

We  now  begin  to  illustrate  how  the  stable  character  of  P„  onset  can  be  used 
to  develop  a  discrimination  algorithm.  The  first  step  involves  measuring  the 
average  waveshape.  To  accomplish  this  we  simply  stack  4  to  6  observations 
for  each  station  and  test  site  This  process  is  illustrated  for  MNV  and  Yucca 
Valiev  events  in  Figure  4.  We  particularly  select  records  which  show  the  split 
svir.g  clearly  sin.'e  this  feature  is  caused  by  pP,  and  we  wish  to  design  our 
d  i  sc  r  i  m  i  t'ar.t  to  detect  it.  It  is  interesting  to  note  the  stability  of  the  feature 
with  respect  to  event  size.  The  largest  event  on  top  is  more  than  2  full  magnitude 
units  larger  than  the  smallest,  yet  the  difference  In  the  character  of  the  "c" 
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Average  Pn  at  MNV 


h  =  640  m 
mK=  5.9 


567  m 


564  m 


564  m 


2  sec 


Sum 


Figure  U.  An  Illustration  of  the  procedure  of  obtaining  the  average 
waveform,  in  this  Instance  at  KNV.  Four  to  six  waveforms  which  exhibit  the 
split  "c"  swing  are  selected,  aligned  in  time  and  summed. 


swing  is  minimal.  The  source  scaling  effect  is  obvious  in  that  the  smallest 
event  is  shorter  in  period  than  the  largest.  It  is  also  worth  noting  that  all 
events  shown  are  deep,  so  the  effect  of  burial  depth  on  pP  character  is  not 
apparent.  We  have  developed  average  P„  waveforms  like  the  one  at  the  bottom  of 
Figure  4  for  all  stations  and  test  sites  where  it  seemed  appropriate,  but  we 
will  not  present  the  details  of  the  averaging  on  a  station  by  station  basis. 
We  will  summarize  the  results  at  the  end  of  this  section. 

KN’B  .Ml  of  the  data  collected  and  analyzed  so  far  are  shown  in  seismic  sections 
in  Appendix  1.  We  will  discuss  these  sections  qualitatively  for  the  next  several 
stations  to  make  a  few  key  points.  The  first  station  will  be  Kanab,  Utah.  At 
this  site  the  structure  of  the  P„  onset  for  Pahute  events  strongly  resembles  that 
obtained  in  Figure  4  by  averaging  the  HNV  signals.  The  Pg  signals  are,  on  the 
other  hand,  completely  different  between  the  two  stations.  At  KNB,  the  Pg appears 
to  be  dominated  by  a  single  pulse.  At  MNV  there  appeared  to  be  several  pulses 
for  Pahute  sources  (see  Figure  2).  Tlie  physics  of  Pg  propagation  has  changed  in 
a  major  way  in  the  two  different  azimuths.  It  will  therefore  be  difficult  to 
develop  a  model  to  explain  how  a  discriminant  based  on  Pg  is  working,  and  its 
transportability  will  be  difficult  to  establish.  Again  there  are  no  other  phases 
visible  in  the  section  except  P„  and  .  Note  that  the  ratio  of  the  two  phases 
is  again  about  5  to  1 . 

The  instability  in  Pg  between  signals  from  Yucca  Flat  and  Pahute  is  even 
more  remarkable  than  that  between  HN\’  and  KNB.  There  is  a  great  increase  in  the 
complexity  of  Pg  though  the  path  differences  for  the  two  test  sites  are  small. 
The  most  logical  explanation  Is  chat  strong  scattering  near  the  source  is  causing 
this  difference.  It  is  difficult  to  understand  how  such  a  strong  effect  would 
not  i  rC  ei  f  ere  with  spectral  di  sci  iminat  ion  us  ing  Pg  .  The  average  P^  onset  ma  intain.s 
i  ' V  ihle  form  and  the  ratio  of  Pj..  to  Pg  remains  the  same. 
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JAS  In  our  earlier  report,  where  we  first  described  and  characterized  the  onset 
of  P„ ,  we  worked  with  data  froii!  the  DWWS5N  stations  Albuquerque  and  Jamestown. 
We  showed  that  the  split  third  swing  was  a  very  common  observation  at  JAS  and 
that  the  average  waveform  was  very  much  like  the  one  shown  in  Figure  4.  The 
scattering  of  Pg  is  very  strong  for  both  test  sites  and  the  amplitude  ratio  of 
the  two  phases  is  5  to  1 . 

ALQ  The  remainder  of  this  section  deals  with  stations  at  which  interpretation 
of  Pn  waveforms  is  somewhat  more  problematical.  At  ALQ  and  ELX,  the  receiver 
structure  seems  to  complicate  the  waveform  slightly.  At  ALQ,  the  first  swing 
of  the  waveform  is  consistently  much  more  emergent  than  at  other  stations.  Other 
than  that,  the  waveshape  is  fairly  consistent  with  those  at  other  stations  in 
the  digital  network.  Figure  5  shows  six  observed  signals  from  Yucca  Flat  events 
with  varying  depth.  The  feature  associated  with  pP  is  indicated  by  arrows.  Its 
timing  with  respect  to  the  first  energy  is  comparable  to  that  observed  at  other 
.stations,  but  the  overtoil  waveform  is  somewhat  different.  The  depth  dependence 
of  the  splitting  feature  is  clear  supporting  the  premise  that  it  is  related  to 
the  free  s^irface  int«.  rac  t  ion .  The  level  of  scattering  of  Pg  is  different  between 
Yucca  F''2t  and  Pahure  Me.sa  again  being  stronger  for  the  former  site.  The  amplitude 
ratio  is  comparable  to  that  at  other  stations. 

ELK  The  situation  at  Elko,  Nevada  is  much  like  that  at  Albuquerque.  The  receiver 
function  is  somewhat  complicated,  but  the  physics  of  the  wave  propagation  seems 
consistent  with  other  site.s.  Figure  6  shows  a  representative  set  of  waveforms 
with  arrows  indicating  a  feature  that  might  be  associated  with  pP.  This  type 
of  v/iriaritui  of  waveform  due  to  receiver  structure  is  not  unexpected  given  that 
t 'n  i  s  i sMcii  hif.h  frequencv  data.  Al  teleseismic  (30°  -  90°)  ranges,  short  period 
'wU'.S'S'.V  p  wa-.o.s  hive  a  reLjri\v>;v  ..on-'dan'  .shape  except  at  a  small  number  of 
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YUCCA  to  ALQ  or  ANMO 
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Figure  5.  Observed  Pj^  signals  from  Yucca  Valley  events  at  DWWSSN  station  ALQ. 
The  split  "c"  swing  is  interpreted  as  the  effect  of  pP^.  As  events  get 
larger,  the  increase  in  the  duration  of  the  source  function  washes  out  the  pP 
effect . 
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Figure  6.  Observed  signals  at  ELK.  The  waveform  shows  additional 
complexity  at  the  beginning.  A  feature  which  seems  event  depth  dependent  and 
might  be  associated  with  pP^,  is  indicated  by  arrows. 


stations  with  complex  receiver  structures.  Such  variations  can  be  removed  by 
relative  waveform  analysis  techniques  (Burger  et  al .  ,  1987).  The  other  aspects 
of  the  regional  P  wave  train  at  ELK  remain  constant. 

PAS  The  reason  that  special  note  has  been  made  in  the  preceding  discussion  of 
the  ratio  of  Pj,  to  Pg  is  that  it  varies  markedly  at  stations  to  the  south  of  NTS. 
The  map  in  Figure  1  shows  that  these  are  at  Landers,  Pasadena  and  Pinyon  Flat 
California.  The  latter  two  have  only  recently  been  upgraded  to  broad  band  digital 
stations  with  state-of-the-art  Streckeisen  instruments.  Because  they  are  rel¬ 
atively  new,  the  available  data  from  them  are  limited.  At  PAS,  the  waveform  of 
P,  onset  is  not  greatly  different  than  at  most  stations  in  the  digital  network. 
Figure  7  shows  the  averaging  of  three  waveforms  to  sum  to  a  final  result  much 
like  that  in  Figure  2  for  MNV.  However,  there  is  a  marked  change  in  the  relative 
amplitude  of  at  PAS  as  shown  in  Figure  8.  On  the  top,  are  shown  three  regional 
P  signals  from  Yucca  events  at  JAS.  It  i.s  illustrated  in  the  appendix  that  these 
are  signals  tN’pical  of  those  from  every  station  discussed  to  this  point  in  terms 
of  the  Pj,  to  Pg  ratio.  The  bottom  three  traces  are  typical  of  those  observed  at 
P.AS .  The  ratio  has  decreased  from  5  to  1  to  less  than  2  to  1.  The  seismic 
sections  in  the  appendix  show  clearly  that  PAS  is  different  in  this  regard  for 
both  Pahute  and  Yucca  Flat  events.  The  reason  is  that  the  Pj,  is  being  amplified 
by  a  rapid  velocity  increase  or  discontinuity  in  the  lid. 

LAC  The  station  at  PAS  is  intermediate  in  range  between  LAC  and  PFO  and  is 
apparently  located  in  the  middle  of  the  triplication.  LAC  is  at  the  beginning. 
The  events  closest  to  LAC  at  a  range  of  about  295  km  have  a  normal  P^.  Pj, 
waveforms  from  several  events  are  compared  at  stations  MNV  and  LAC  in  Figure  9. 
The  correspondence  between  waveforms  is  very  strong  indicating  a  simple  and 
deterministic  type  of  wave  propagation.  The  scaling  of  frequency  content  with 
I  V.  Ti*  .size  is  very  consistent  between  the  two  stations.  At  slightly  greater 
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Figure  7.  The  average  waveform  at  PAS.  The  "c"  swing  is  consistently 
split  even  though  the  P^  is  being  amplified  due  to  a  triplication. 


Figure  8.  Amplification  of  with  respect  to  P  at  PAS.  At  JAS  and  all  other 
stations  not  to  the  south  of  NTS,  the  ratio  of  the  phases  is  about  5  to  1 .  At 
PAS  and  other  southern  stations  it  is  less  than  2  to  1 . 


LAC 


MNV 


Figure  9.  A  comparison  of  waveforms  for  common  events  at  MNV  and  LAC.  At 
these  ranges  the  triplication  has  not  begun  to  have  an  effect  and  the 
waveshapes  are  very  similar. 


ranges,  the  triplication  clearly  emerges,  however,  as  shown  in  Figure  10.  By 
about  315  km,  a  clear  precursor  has  emerged  ahead  of  a  strong  second  arrival. 
The  section  in  the  appendix  shows  that  grows  systematically  with  respect  to 
Pg  though  not  quite  to  the  levels  it  reaches  at  PAS  farther  into  the  triplication. 
PFO  As  shown  in  the  final  sections  in  the  appendix,  the  data  available  from 
Pinyon  flat  is  very  limited.  The  waveform  of  the  is  relatively  complicated, 
and  the  P,,  to  Pg  ratio  decreases  with  range.  We  interpret  this  as  the  end  of  the 
triplication.  PFO  may  also  have  a  complex  receiver  function,  but  resolving 
turning  point  from  near  receiver  structural  effects  is  beyond  the  scope  of  this 
investigation. 

To  quantify  the  amplification  effect,  we  have  measured  the  to  Pg 
amplitude  ratio  on  a  set  of  the  available  records  from  the  California  stations. 
The  P^  arr.plitude  is  defined  as  the  RMS  value  within  5  s  after  the  first  arrival. 
The  Pg  amplitude  is  defined  as  the  RMS  value  within  20  s  after  the  arrival  time 
of  PmP  as  predicted  by  the  Helmberger  and  Engen  (1980)  crustal  model.  This  model 
appears  from  experience  to  be  a  good  predictor  of  Pg  onset.  Since  it  has  been 
made  clear  that  the  behavior  of  Pg  varies  between  the  Pahute  Mesa  and  Yucca  Flat 
test  sites,  it  would  be  best  to  deal  with  data  from  only  one  site.  More  data 
are  available  from  the  Yucca  Flat  site,  so  at  LAC,  and  PAS  we  have  used  data 
from  only  that  site.  Very  little  data  are  available  from  either  site  at  PFO. 
Therefore,  we  used  data  from  both  test  sites  at  this  station.  The  resulting 
amplitude  ratios  are  plotted  as  a  function  of  distance  in  Figure  11.  Certain 
features  of  the  behavior  are  very  clear,  such  as  the  relative  amplification  of 
Pj.,  at  PAS  and  PFO.  The  upward  trend  at  LAC  and  downward  trend  at  PFO  are  more 
questionable  but  are  suggested  by  dashed  lines  in  the  figure. 
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Figure  10.  Development  of  the  triplication  with  range  at  LAC.  At  the  close 
ranges,  the  first  energy  is  the  largest  and  the  waveform  is  similar  to  that 
at  most  stations.  At  greater  ranges,  a  large  second  arrival  develops  and 
moves  toward  the  front  of  the  record. 


Ratio 


o 


Figure  11.  The  ratio  of  the  RMS  aniplltude/sec ,  of  to  with  range.  It  is 
very  clear  that  the  ratio  increases  between  LAC  and  PAS.  Other  trends  which 
rr.ight  be  significant  are  indicated  by  dashed  lines. 
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Results  of  the  Pj,  Averaging  Despite  the  preceding  discussion  of  some  of  the 
complexities  of  P,,  propagation  in  the  western  U.  S.  ,  we  believe  that  the  importance 
of  the  phase  is  its  stability  and  simplicity  in  most  areas.  Triplication;  due 
to  structure  in  the  lid  are  not  a  universal  of  the  earth  feature  or  they  would 
have  been  noticed  ’ ang  before  The  existence  of  one  in  this  region  is  probably 
related  to  th<?  intense  tectonics  of  southern  California.  Such  phenomena  are 
less  likely  to  occur  in  stable  continent.  We  expect  that  there  will  always  be 
some  stations  with  complex  receiver  structures  such  as  ELK,  but  if  the  dis¬ 
crimination  procedure  is  based  on  obsem’at ions  from  an  array  of  stations  this 
di f { icul tv  can  be  overcome . 

Ti'.e  tes's  at  Pcihv;te  .Mesa  are  on  average  larger  than  those  from  Yucca  Flat. 
T't'e  striK.tuie  there  alse  ap-.u.c.r,s  to  b<'  simpler  and  the  scattering  less  intense, 
rh-.  ■,  t'  efii  .'t  .Sei  ms  in  genet. i'.  tt  appear  more  consistently  and  more  clearly  for 
!  ivtuis.  r'igure  11’  sh  'v.;  the  average  P^  for  fi  stations  in  the  digital 

r.e'wor’r'.  The  splitting  t  iit'  ’c"  svir^g  is  indicated  by  arrows.  The  only 
stations  in  the  digital  n.etworK  tor  n-hich  ve  do  not  show  a  result  are  E!J< ,  LAC 
at'.d  FFO  I'he  former  h.as  ti  coi.uif't  re. reiver  strec-tur*',  and  the  latter  two  are 
aft.  te  i  hv  'tie  rri;  licaticn.  Figure  1  i  ‘.'rows  similar  averaging  results  for  7 
St  i'  :  ru;  .  Ail  of  t!ie  di  it.L.  st/r.;.  .::c.,-p:  FiLK  and  PFO  are  repre.ser.ted  Ve 
i.  -htri'!.  t  Iv  ;<-su!'  f  I.'-  tiv  rag  i  tij  sigiials  at  ranges  prior  to  the 

’  ■  :  ,  ni . :  t!; ,  c  ;  .  ; f  r  i  ;  ■  1  i  a  t  ;  ■  n  V.-  .-is-  -.h.-v  j .  .:n  !  t  r  f  i  .  .r  two  an.i  1  og  stations, 
r.r''.  csi'  A  .  Wi  d’t  t  in  e  ;;  pt.  ■if.  :.,  report  .  we  hand  digitized  the  data 

i  1 'c;  ‘  ■  ,i :  i  cr-s  .rd  t  !.  :  ■ '■,  t  hr  s.o:'.-  w.iv  C’ven  the  stabiiitv  of 

*T'  -..i"-  ‘  i  :■  '  i,.  t  „■  t’r  it*  .utf  w.-jjd  V:.r-e  to  roiM  lude  tha’  w  c  are  seeing 

!  r  ;  f  e '  ■  •  1  or,  (i‘  i  :  tn-  \]  '.rop-rtv  >f  explosion  sources.  V'e  will 

■'.■■'W  in  •  ic'  tc'  inw’.;;,',  'h-'  t  >  •.  .  ;  .o;,,,;  ir:Vf,l  '(■  .  :  t- 1  1  e  c  '  1  r.  n  of  e’^ergv  from  these 
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Figure  12.  The  average  waveform  for  Pahute  events  observed  at  stations  in 
tVie  western  U.S.  digital  net.  The  split  third  swing  is  indicated  by  arrows. 


Figure  13.  The  average  waveforrn  for  Yucca  events  obser\’ed  at  stations  1 
the  western  U.S.  digital  net.  The  split  third  swing  is  indicated  by  arrows. 


MODELING  OF  THE  WAVEFORMS 


In  the  preceding  discussion,  we  have  dealt  with  signals  with  a  uniform  WWSSN 
short  period  response.  In  the  modeling  studies  to  be  presented  next,  we  will 
work  with  two  of  the  LLNL  stations  where  the  pP  effect  is  very  clear.  We  will 
take  advantage  of  the  very  broad  band  response  of  the  LLNL  instruments  and 
deconvolve  the  instrument  response  cut.  We  will  work  with  actual  ground  velocity 
traces  and  show  that  the  effect  of  pP  is  simple  and  straightforward  to  understand. 

Tfie  theoretical  instrument  response  function  for  the  LLNL  stations  was 
developed  and  presented  by  Vergino  et  al .  (1986).  They  also  discussed  the 
procedure  for  deconvolving  the  response  out  to  generate  true  ground  velocity. 
Cln  this  instance,  deconvolution  simply  implies  inverse  convolution  whereas  it 
has  been  used  in  other  cases  to  imply  more  complex  procedures.)  With  the 
instrument  effect  removed,  we  can  represent  an  observed  velocity  trace  V(t)  in 
the  standard  synthetic  body  wave  decomposition  as 

V(t)  -  S(t)  *  G(t)  *  Q(t) 

wViere  S  is  the  source,  G  is  the  elastic  response  of  the  earth  and  Q  is  a  correction 
for  ane last ic i ty .  The  *  represents  convolution.  For  the  purposes  of  this 
discussiori  it  is  convenient  to  further  decompose  the  elastic  response  G(t)  into 

G(t)  -  G„(t)  *  Gt(t)  *  G,(t) 

These  three  subresponses  represent  the  effect  of  structure  near  the  source, 
turning  point  and  receiver  respectively.  Such  a  decomposition  is  justified  for 
iViorr  '.imes  near  the  first  arrival.  Gg(t)  contains  the  effects  of  the  free 
surface  and  near  source  crust.  Gj(t)  contains  the  effects  of  structure  near  the 
rurnir.g  point  if  any  or  the  effect  of  the  ray  turning  process  Itself.  At  complex 
receiver  .sites  such  as  EIJ< ,  Gj.(C)  would  have  a  major  effect.  In  the  following, 
wf  siirll  mode’  data  from  the  simple  sites  at  MNV  and  KNB.  Ve  will  neglect  the 
f  if.  •  ■.  cif  tin's  la.st  operator.  Each  of  the  operator.s  in  this  general  decomposition 
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have  been  studied  and  characterized  extensively  in  the  past  but  apparently  not 
using  the  regional  P„  waveform.  The  data  set  assembled  in  this  study  presents 
a  unique  opportunity  to  do  so. 

The  most  recent  work  on  the  nuclear  explosion  source  function  S(t)  by  Murphy 
(1988)  indicates  that  it  can  be  well  approximated  by  the  Mueller  Murphy  repre¬ 
sentation  for  Pahute  wet  tuff  for  all  test  sites  providing  the  events  are  well 
coupled  and  below  the  water  table.  There  is  a  need  for  further  testing,  however, 
of  both  whether  the  functional  form  of  the  Mueller  Murphy  source  is  appropriate 
and  whether  the  yield  scaling  is  supported  by  the  data.  The  common  assumption 
regarding  Gt^(t)  is  that  it  can  be  represented  at  short  times  by  a  step  function. 
This  corresponds  to  the  standard  wave  front  expansion  from  ray  theory  for  a  head 
wave  (Aki  and  Richards,  1980).  In  the  following,  we  will  explore  the  possibility 
that  it  might  be  better  represented  by  a  delta  function.  In  terms  of  the  same 
wave  front  expansion  theory,  this  would  correspond  to  the  P^  onset  actually  being 
energy  turning  in  the  lid.  That  is  to  say  that  rather  than  velocity  being 
constant  beneath  the  Moho  it  increases  gradually  with  depth.  Forward  calculations 
show  that  even  a  weak  gradient  will  cause  the  turning  point  response  to  evolve 
into  a  delta  function  rather  than  a  step  function  very  quickly.  The  question 
of  how  strong  the  correction  for  attenuation  should  be  has  been  debated  for  many 
years  for  both  teleseismic  and  regional  P  waves  (Burger  et  al . ,  1987).  It  has 
generally  been  observed  that  both  and  follow  a  relatively  high  Q  path. 

In  the  following,  we  will  present  accurately  computed  synthetics  for  realistic 
crystal  structures  and  compare  them  to  observations.  However,  we  begin  with  a 
si'rple  .schejuat  ic  calculation  to  illustrate  the  features  of  most  importance  in 
the  data  aiui  the  synthetics.  Figure  14  shows  this  calculation.  The  top  trace 
is  the  Mueller  Murphy  source  function  for  velocity  assuming  a  turning  ray.  It 
hat,  been  convol  /ed  with  a  Futterman  operator  with  a  t*  level  set  at  0.1  s.  If 
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Figure  14.  A  simplified  synthetic  calculation  indicating  the  expected  effect 
of  elastic  pP^  on  the  ground  velocity  trace.  It  lengthens  the  downswing  and 
adds  a  second  upswing.  An  observed  trace  from  HNV  Is  shown  at  the  bottom. 
The  features  associated  with  P  and  pP  in  the  data  are  Indicated  by  arrows. 


a  head  wave  model  for  is  assumed  this  pulse  would  simply  be  integrated.  The 
second  trace  represents  elastic  pPn-  It  Is  reversed  in  sign,  delayed  and  reduced 
in  amplitude.  In  this  simplified  calculation,  the  resulting  synthetic  in  the 
chira  row  is  just  the  sum  of  the  first  two  traces.  The  direct  P„  results  in 
first  an  upswing  and  then  a  downswing.  The  pP^  results  in  the  opposite.  The 
effect  of  pPn  is  to  cause  a  shoulder  in  the  P^  and  to  extend  its  duration.  The 
observed  P^  ground  velocity  from  Yucca  flat  event  JORNADA  is  shown  at  the  bottom. 
The  features  associated  with  the  direct  and  pP  arrivals  are  indicated  by  arrows. 

The  calculation  shown  in  Figure  14  is  highly  approximate  in  that  the  delay 
time  and  relative  amplitude  of  pP^  have  been  set  arbitrarily.  These  values  along 
with  the  effects  of  crustal  structure  should  be  computed  from  an  accurate  model 
of  the  NTS  crust.  In  the  following,  we  compute  the  response  Gg(t)  using  standard 
layer  matrices  and  an  assumption  of  constant  slowness.  The  crust  model  we  use 
is  based  on  modeling  studies  of  near  field  data  and  has  been  shown  to  be  accurate 
for  timing  purposes  (Hartzell  et  al.,  1983).  We  use  the  appropriately  scaled 
Mueller  Murphy  source  for  S(t)  and  again  a  Futterman  operator  for  Q(t). 

Figures  15  and  16  compare  data  to  synthetics  for  four  Pahute  Mesa  events 
as  observed  at  MNV .  From  these  calculations,  we  hope  to  address  a  number  of 
questions.  These  include,  is  the  turning  ray  or  head  wave  model  for  G^Ct)  more 
appropriate,  does  the  Mueller  Murphy  source  scale  correctly  and  what  is  an 
appropriate  value  of  t*  for  ?„?  We  also  would  like  to  know  whether  the  commonly 
used  approximation  of  pi  as  an  elastic  reflection  of  P  is  accurate. 

The  results  for  the  turning  ray  approximation  of  G^Ct)  are  shown  in  Figure 
15  and  for  the  head  wave  in  16.  We  began  by  simply  assuming  that  the  Mueller 
Mill  tihy  soun  e  gives  the  correct  source  pulse  shape  S(t).  Because  the  Mueller 
■Murphy  .source  assumes  an  f'-^  spectral  decay  rate,  the  predicted  velocity  pulse 
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Figure  15.  Observed  versus  synthetic  wavefori 
The  turning  ray  model  for  propagation  is  assumed 
scaling  relations  and  a  t*  value  of  0.1  s  are  al 
arrows,  the  observed  pP  arrival  is  late  but  appro 
the  elastic  prediction. 
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Figure  16.  Observed  versus  synthetic  waveforms  at  MNV  for  Pahute  events.  All 
assumptions  are  the  same  as  for  Figure  15  except  for  a  head  wave  model  for 
propagation. 


is  discontinuous.  We  increased  the  level  of  attenuation  (t*)  until  the  dis¬ 
continuous  pulse  was  smoothed  into  approximately  the  observed  pulse.  There  are 
two  features  of  the  observed  pulses  to  be  considered.  The  first  is  the  zero  to 
peak  or  rise  time.  Ine  second  is  the  duration  of  the  first  swing.  The  t*  value 
for  the  synthetics  in  Figures  15  and  16  was  ultimately  set  at  0.2  s.  The 
comparisons  in  Figure  15  show  that  for  the  turning  ray  model,  this  value  gives 
approximately  the  right  duration  for  the  first  swing  but  that  the  theoretical 
rise  tim.e  remains  too  high.  The  latter  discrepancy  could  mean  that  the  f"^ 
spectral  decay  assumption  is  incorrect  or  that  a  frequency  dependent  Q  correction 
is  required.  The  Mueller  Murphy  source  scaling  predicts  a  larger  difference  in 
the  durations  of  the  largest  and  smallest  events  than  is  observed.  This  is  a 
small  data  sample,  but  it  would  appear  that  the  possibility  that  the  scaling 
effect  is  too  strong  in  the  Mueller  Murphy  formalism  will  need  further  exploration 
in  the  future.  The  computed  pP  arrivals  and  what  we  interpret  as  pP  arrivals 
in  the  data  are  indicated  by  arrows.  The  observed  pP  is  clearly  and  consistently 
later  than  the  predicted.  This  phenomenon  of  anomalously  delayed  pP  has  been 
observed  consistently  in  teleseismic  data  in  many  different  investigations 
involving  many  different  approaches  (Murphy,  1988;  Lay,  1985).  The  P^  waveforms 
confirm  this  anomaly.  The  only  interpretation  warranted  is  that  what  appears 
as  pPj,  or  teleseism.ic  pP  is  more  complex  than  a  simple  elastic  reflection.  The 
observed  amplitude  and  waveshape  of  pP  in  the  figure  is  comparable  to  the 
predictions.  Unfortunately,  the  traces  had  to  be  truncated  after  only  two  seconds 
because  of  the  onset  of  Pg  (see  Figure  2). 

The  head  wave  synthetics  in  Figure  16  do  not  fit  the  observations  nearly 
as  well  as  the  turning  ray  synthetics.  Because  of  the  additional  integration, 
t}if  svnrhetic  pulses  become  much  too  long  in  period.  The  shoulder  associated 
with  pP„  onset  washes  out  completely.  The  value  of  t*  was  adjusted  to  match  the 
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observed  rise  time  of  the  pulses.  As  the  figure  indicates,  it  is  impossible  to 
simultaneously  match  the  rise  time  and  the  duration  of  the  first  peak  while  using 
the  head  wave  model  for 

Similar  comparisons  between  data  and  synthetics  for  Yucca  Flat  events 
observed  at  MNV  are  shown  in  Figures  17  and  18.  In  this  instance,  a  t*  value  of 
0.1s  appears  most  consistent  with  the  data.  The  match  between  observations  and 
synthetics  in  the  top  two  traces  of  Figure  17  is  extraordinary.  There  can  be 
little  doubt  that  a  feature  associated  with  the  free  surface  interaction  is  being 
modeled.  It  is  again  clear,  however,  that  this  feature  arrives  too  late  to  be 
elastic  pP^.  It  also  appears  that  in  some  cases  pP,,  has  a  larger  amplitude  than 
direct  P^.  This  could  only  be  caused  by  some  nonlinear  effect  in  the  source 
region.  Under  the  circumstances ,  it  is  probably  best  to  refer  to  the  pulse  as 
effective  pP^  and  not  to  attempt  to  model  it  too  closely  using  assumptions  of 
perfect  elasticity.  The  pP^,  from  the  magnitude  5.^  event  at  the  lower  left  shows 
additional  structure.  This  is  another  indication  that  a  complex  process  is 
occurring . 

The  Mueller  Murphy  source  scaling  again  fails  to  predict  variations  in 
frequency  content  associated  with  source  size.  There  is  a  large  shift  in  frequency 
content  between  the  top  events  with  magnitude  5.7  and  5.9  which  the  scaling  law 
does  not  predict.  There  is  a  modest  increase  in  frequency  content  as  source 
size  decreases  from  mj,  5.7  to  5.4  to  4.5,  but  the  scaling  laws  over  predict  the 
effect.  The  head  wave  model  synthetics  in  Figure  18  are  uniformly  too  low  in 
frequency  content.  This  reinforces  the  credibility  of  the  turning  ray  repre¬ 
sentation  of  the  Pn  arrival. 

Similar  comparisons  of  data  from  Pahute  and  Yucca  Valley  at  station  KNB  are 
shown  in  Figures  19,  20,  21  and  22  respectively  for  the  same  events  modeled  at 
MNV.  We  will  not  discuss  these  comparisons  in  extended  detail,  but  we  believe 
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Figure  17.  Observed  versus  synthetic  waveforms  at  MNV  for  Yucca  events.  All 
assumptions  are  the  same  as  for  Figure  15  except  for  an  increase  in  t*  from 
0.1  to  0.2  s . 


Figure  18.  Observed  versus  synthetic  waveforms  at  MNV  for  Yucca  events.  All 
assumptions  are  the  same  as  for  Figure  15  except  for  a  head  wave  model  for 
propagation  and  an  increase  in  t*  from  0.1  to  0.2  s. 
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Figure  ?1.  Obser\'ed  versus  synthetic  waveforms  at  KNB  for  Yucca  events, 
assumptions  are  the  same  as  in  Figure  17. 
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thev  illustrate  the  following  points.  Effective  pP„  is  always  present  as  indicated 
by  the  arrows.  It  is  always  significantly  later  than  the  elastic  predictions. 
For  Pahute  events,  the  amplitude  of  effective  pP  is  approximately  the  same  as 
for  the  elastic  predictions.  For  Yucca  Valley  events  it  is  consistently  larger. 
The  physics  of  effective  pP  generation  appears  to  be  test  site  dependent  The 
value  of  t*  appears  to  be  about  0.1  or  0.2  s  for  P^j.  This  is  almost  an  order  of 
magnitude  decrease  from  the  value  for  teleseismic  P.  The  head  wave  model  for 
G,(t)  consistently  produces  synthetics  which  are  too  long  period  to  fit  the 
obser'.’ations .  The  Mueller  Murphy  scaling  laws  do  not  successfully  predict  the 
observed  frequencv  shifts  for  Yucca  Valley  events. 
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DISCRIMINATION  WITH  THE  WAVEFORM 

Figures  12  and  13  demonstrate  that  the  waveform  is  on  average  very  stable 
and  that  it  exhioics  a  distinctive  feature  a  few  seconds  after  onset.  The  above 
modeling  studies  provide  a  convincing  case  that  this  feature  is  associated  with 
effective  pP^.  Since  earthquakes  uniformly  occur  much  deeper  than  explosions, 
it  is  reasonable  to  expect  that  few  earthquakes  will  exhibit  a  highly  similar 
feature  at  such  short  times  into  the  waveform.  In  a  previous  report  (Burdick 
et  al . ,  1988),  we  demonstrated  in  a  quantitative  way  that  this  was  the  case. 
We  collected  P„  waveforms  from  the  earthquakes  mapped  in  Figure  1  at  stations 
ALQ  and  JAS .  We  windowed  out  3  second  onset  waveforms  from  both  the  earthquake 
and  explosion  signals  and  measured  their  correlation  with  an  analytic  norm.  This 
norm  has  been  described  in  detail  by  Burdick  (1986)  .  It  attains  a  maximum  value 
of  1.0  if  and  only  if  two  waveforms  are  identical.  Figure  23  shows  the  result 
of  measuring  the  correlation  of  Pahute  explosion  versus  earthquake  waveforms 
with  the  average  explo.sion  waveform  from  ALQ-  Pahute  signals  seem  to  be  more 
consistent  in  their  pPj,  character  than  Yucca  signals.  The  discrimination  of 
explosion  from  earthquake  signals  is  very  good  at  all  stations  just  as  in  Figure 
23.  Unfortunately,  the  majority  of  events  detonated  at  Pahute  Mesa  are  large. 
Standard  techniques  could  be  used  to  discriminate  them  from  earthquakes.  Smaller 
magnitude  explosions  are  usually  detonated  at  Yucca  Flat,  and  because  of  the 
greater  instability  of  pP^,  discrimination  is  somewhat  more  problematical  there. 

Figure  24  shows  the  results  for  Yucca  Flat  events  at  ALQ.  The  discrimination 
performance  of  the  waveform  norm  definitely  deteriorates  between  magnitude  4.0 
and  4.5.  However,  a  careful  review  of  this  deterioration  shows  that  it  is  caused 
at  least  in  part  by  a  signal  to  noise  ratio  problem.  ALQ  is  the  most  distant 
St, it  ion  in  the  digital  net,  so  this  might  not  be  so  much  of  a  problem  at  other 
.-t.ntions.  Also,  it  generally  seems  that  P„  signal  to  noise  ratio  is  much  better 
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Figure  23.  Discrimination  of  Pahute  explosions  from  earthquakes  using 
correlation  with  the  average  waveform  at  ALQ.  The  explosions  are  stars  and 
the  earthquakes  squares. 
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Figure  24.  Discrimination  of  Yucca  explosions  from  earthquakes  using 
correlation  with  the  average  wavefora  at  ALQ.  The  explosions  are  stars  and 
the  earthquakes  crosses. 
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in  stable  continental  regions  than  in  the  highly  tectonized  western  U.S.  Figure 
2  5  shows  the  discrimination  results  for  Yucca  Flat  events  at  JAS .  The  separation 
of  explosions  from  earthquakes  between  4.0  and  4.5  is  somewhat  improved. 

A  further  consideration  in  developing  a  discriminant  based  on  waveform 
is  that  some  dependence  of  the  waveform  on  event  size  should  be  expected.  In 
the  examples  of  the  waveform  discriminant  shown  so  far,  we  have  simply  tested 
for  the  correlation  of  waveforms  with  the  average  waveforms  shown  in  Figures  12 
and  13.  Generally  those  averages  were  formed  from  signals  of  larger  events  since 
thev  tend  to  exhibit  the  splitting  of  the  third  swing  more  clearly.  To  estimate 
the  importance  of  event  size  on  the  waveform  discriminant ,  we  averaged  explosion 
P,,  waveforms  at  JAS  from  six  explosions  with  magnitude  between  4.0  and  5.0.  The 
correlation  of  that  waveform  with  the  JAS  discrimination  data  base  is  shown  in 
Figure  26.  The  correlation  of  the  large  explosions  with  the  average  waveform 
elea,.’y  with  respect  to  the  levels  Figure  25.  The  discrimination  efficiency 
at  small  event  size  clearly  improves.  It  appears  that  as  the  waveform  discriminant 
is  developed  further,  an  event  size  correction  should  be  incorporated  into  it. 

The  discrimination  tests  showT.  so  far  demonstrate  only  that  explosion 
waveforms  are  consistent  and  stable  at  regional  stations  and  systematically 
different  than  earthquakes.  The  conditions  of  the  tests  are,  however,  highly 
idealized.  In  a  realistic  situation,  signals  will  be  recorded  on  a  regional 
array  from  continuously  varying  azimuths.  The  discrimination  procedure  would 
have  to  assimilate  all  of  the  waveforms  from  the  array  and  test  whether  the 
signals  had  a  feature  resembling  explosion  pP^  on  average.  The  consistency  of 
the  Pf  wa'.^eforms  in  Figures  12  and  13  indicates  that  such  a  procedure  should  be 
effective  at  some  level.  A  .simple  test  as  to  whether  explosion  pP^  signals  from 
or.,  st.ition  can  be  used  to  discriminate  between  events  at  another  station  has 
b'-' '  r  nU'  by  coMelating  the  average  MNV  .signal  with  the  JAS  data  set.  The 
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Figure  25.  Discrimination  of  Yucca  explosions  from  earthquakes  using 
correlation  with  the  average  P„  waveform  at  JAS .  The  explosions  are  stars  and 
the  earthquakes  are  crosses. 
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Figure  26.  Discrimination  of  Yucca  explosions  from  earthquakes  using 
correlation  with  an  average  P„  waveform  explicitly  formed  by  averaging  signals 
from  smaller  events.  In  comparison  to  Figure  25,  discrimination  is  better  for 
small  events  and  poorer  for  large  events.  This  is  presumably  an  effect  of 
source  scaling. 
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results  are  shown  in  Figure  27.  The  discrimination  efficiency  is  comparable  to 
the  level  achieved  with  the  JAS  waveform  proving  the  feasibility  of  the  waveform 
discriminant . 
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Figure  27.  Discrimination  of  Yucca  explosions  from  earthquakes  using  the 
average  waveform  from  MNV  and  the  data  base  from  JAS.  Explosions  are  stars 
and  earthquakes  are  crosses.  The  separation  of  the  populations  is  comparable 
to  that  in  Figure  26  illustrating  the  transportability  of  the  discriminant. 
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DISCUSSION 


ihere  are  balancing  concerns  involved  in  the  development  of  an  effective 
regional  discrimination  procedure.  On  the  one  hand,  a  high  degree  of  reliability 
is  desired,  but  on  the  other  so  is  a  dependable  level  of  transportability.  In 
the  western  U.S. ,  the  spectral  discriminant  of  Taylor  et  al.  (1988)  seems  to  be 
more  effective  than  the  waveform  discriminant  as  it  has  been  developed  here. 
The  difficulty  with  the  former  is  that  the  physics  behind  it  is  not  understood. 
Most  lines  of  evidence  involving,  for  instance,  near  field  or  teleseismic  data 
indicate  that  there  is  no  consistent  spectral  difference  between  explosions  and 
earthquakes.  Thus  there  is  no  way  to  be  sure  that  the  processes  which  cause 
this  difference  in  the  western  U.S.  will  occur  consistently  in  other  regions. 
The  Pj,  waveform  discriminant  is  based  on  detecting  the  effects  of  the  interaction 
of  seismic  energy  with  the  free  surface.  Albeit  these  interactions  are  apparently 
nonlinear  and  not  completely  understood,  they  will  occur  for  any  nuclear  test. 
The  Pn  waveform  discriminant  can  be  reliably  transported. 

Further  development  of  the  waveform  discriminant  should  involve  the  col¬ 
lection  and  analysis  of  performance  statistics  in  the  western  U.S.  How  often 
does  a  signal  from  a  western  U.S.  earthquake  correlate  highly  with  a  Pj,  signal 
from  an  explosion?  It  could  be  expected  that  similar  performance  statistics 
w'ould  hold  in  a  region  where  actual  verification  was  taking  place .  Some  adjustment 
to  the  discrimination  procedure  and  to  the  performance  statistics  could  be  based 
on  studies  of  P^  waveforms  from  earthquakes  in  more  normal  tectonic  regimes  than 
the  western  U.S.  The  actual  implementation  of  the  discriminant  on  a  computer 
could  follow  one  of  several  approaches.  One  would  be  an  empirical  one  in  which 
observed  waveforms  from  a  regional  net  are  continuously  correlated  with  the 
average  explosion  Pn  waveforms  from  the  western  U .  S  .  digital  net.  Another  approach 
could  be  based  on  synthetic  modeling.  The  models  could  be  adjusted  to  fit  the 


46 


western  U.S.  data  base  and  then  used  to  predict  waveforms  for  other  regions. 
Corrections  for  source  scaling  and  for  variations  in  attenuation  would  be 
straightforward  to  implement  in  this  approach.  A  third  method  could  be  based 
on  recognition  of  patterns  in  the  waveform.  For  example,  the  method  could  be 
based  on  logical  identification  of  the  "c"  swing  and  a  logical  test  for  whether 
splitting  was  occurring. 
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CONCLUSIONS 


The  waveform  of  the  onset  of  has  been  found  to  be  extraordinarily  stable 
in  character  in  the  western  U.S.  The  waveform  for  nuclear  explosions  c ,n  be 
easily  synthesized  using  standard  models  and  methods.  The  specific  modeling 
studies  described  here  demonstrate  that  at  high  frequency  is  more  probably  a 
turning  ray  from  a  gradient  in  the  lid  than  a  true  head  wave  along  the  Moho 
discontinuity.  The  value  of  t*  for  P„  appears  to  be  in  the  range  of  0.1  to  0.2 
seconds.  This  means  that  the  explosion  source  process  can  be  observed  in  much 
closer  detail  than  it  can  from  teleseismic  data  where  t*  values  are  closer  to 
1.0.  The  scaling  laws  for  the  Mueller  Murphy  source  do  not  explain  the  observed 
scaling  of  P„  waveforms  adequately.  The  effective  pP^  can  be  observed  clearly 
in  the  waveform.  It  differs  in  character  between  Pahute  Mesa  and  Yucca  Flat 
event.s.  In  all  cases,  it  appears  that  simple  elastic  reflection  from  the  free 
surface  is  a  poor  model  for  effective  pPj,.  The  arrival  time  is  always  late  with 
respect  to  direct  P^,  and  for  Yucca  Flat  the  pP^  is  larger  than  direct  P,, .  The 
waveform  of  pP„  varies  substantially  in  many  cases  from  the  direct  arrival.  It 
has  been  demonstrated  here  that  the  waveform  of  P^  can  be  used  as  the  basis  of 
a  discrimination  procedure. 
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APPENDIX  Ij. 

SEISMIC  SECTIONS  OF  NTS  EVENTS  FROM  THE  WESTERN  U.S  DIGITAL  ARRAY 
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